The information about the state of a qubit gained by a weakly coupled detector 
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We analyze the information that one can learn about the state of a quantum two-level system, 
i.e. a qubit, when probed weakly by a nearby detector. In particular, we focus on the case when 
the qubit Hamiltonian and the qubit's operator being probed by the detector do not commute. 
Because the qubit's state keeps evolving while being probed and because the measurement data is 
mixed with a detector-related background noise, one might expect the detector to fail in this case. 
We show, however, that under suitable conditions and by proper analysis of the measurement data 
useful information about the state of the qubit can be extracted. It turns out that the measurement 
basis is stochastically determined every time the experiment is repeated. We analyze in detail 
the probability distributions that govern the choice of measurement bases. We also analyze the 
information acquisition rate and show that it is largely unaffected by the apparent conflict between 
the measurement and intrinsic qubit dynamics. We discuss the relation between our analysis and 
the stochastic master equation that describes the evolution of the qubit's state under the influence of 
measurement and decoherence. In particular, we write down a stochastic equation that encompasses 
the usual stochastic master equation for the evolution of the qubit's density matrix and additionally 
contains the measurement information that can be extracted from the observed signal. 



I. INTRODUCTION 

Quantum measurement is a subject of importance both 
as a central part of quantum theory [ij] and for its rele- 
vance to quantum information theory 0. At the most 
basic level, the measurement is commonly thought of as 
a sudden projective operation that leaves the measured 
system in an eigenstate of the probed operator. It should 
be kept in mind, however, that the measurement process 
is a physical process that follows from a physical inter- 
action between the measured system and the measuring 
device. As a result, any realistic measurement is carried 
out over a finite duration. One can therefore have the 
situation where the measurement process occurs concur- 
rently with other physical processes. This is the situation 
that we consider in this paper. 

Although the analysis and results of this paper should 
be valid for a variety of different physical settings, we find 
it useful to formulate the analysis in terms of a specific 
setup. In this paper, we use the example of a double- 
quantum-dot charge qubit that is probed by a quantum 
point contact (QPC) i J, S i 0, i, i, [M [ll [ill, [II, Q 
[isl . [l6t . The principle of the measurement is the fact that 
the current in the QPC is sensitive to the charge state of 
the qubit. By measuring the current passing through the 
QPC, one can infer the state of the qubit. One limitation 
that arises in practical situations is that, in order to mini- 
mize the effects of the detector on the qubit at times when 
no measurement is performed (e.g. during the implemen- 
tation of quantum gates), the qubit-detector coupling is 
set to a value that is small compared to the qubit's en- 
ergy scale. As a result one must deal with some form 
of weak-measurement regime. This type of weak, charge- 
sensitive readout works well when the qubit is biased such 



that the charge states are eigenstates of the Hamiltonian 
and therefore do not mix during the measurement. In 
this case one can allow the detector to probe the qubit 
for as long as is needed to obtain a high signal-to-noise 
ratio, without having to worry about any intrinsic qubit 
dynamics. 

In contrast to the simple situation described above, 
when the detector weakly probes the charge state of the 
qubit while the Hamiltonian induces mixing dynamics be- 
tween charge states, it becomes unclear how to interpret 
a given measurement signal. Given the fact that the sig- 
nal typically contains a large amount of detector-related 
noise and the measurement causes unavoidable relaxation 
in the qubit, it might seem that this type of measurement 
cannot be used to determine the initial state of the qubit, 
i.e. at the time that the experimenter decides to perform 
the measurement. Indeed, there have been a number of 
studies analyzing the measurement-induced decoherence 
and the evolution of the qubit's state in this situation 
d [1, [1, [i, [3, [1, [11, but not the question of how to take 
the measurement data and extract from it information 
about the initial state of the qubit. This question is a 
key issue for qubit-state readout and is the main subject 
of this paper. 

We shall show below that high-fidelity measurement in- 
formation can be extracted from the measurement data, 
provided that additional decoherence mechanisms are 
weak and the readout signal can be monitored at a suf- 
ficiently short timescale. It turns out that not only 
the measurement result, but also the measurement ba- 
sis, is determined stochastically in this case. In spite of 
the uncontroUability of the measurement basis, the mea- 
surement results contain meaningful information about 
the initial state of the qubit. In order to demonstrate 
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this statement, we show how these results can be used 
to perform quantum state tomography on the qubit. 
These results show that under suitable conditions and 
by proper analysis of the measurement data useful in- 
formation about the state of the qubit can be extracted 
from this data. 

This paper is organized as follows: In Sec. II we intro- 
duce the theoretical model for describing a charge qubit 
probed by a QPC. In Sec. Ill we formulate the theo- 
retical framework for the qubit evolution conditioned on 
the QPC's output signal and how measurement informa- 
tion can be extracted from this signal. We also intro- 
duce an equation that supplements the stochastic mas- 
ter equation for purposes of interpreting the measured 
QPC signal. In Sec. IV we present results of numerical 
calculations that display a number of interesting results 
in this system, most importantly the stochastic nature 
of the measurement basis and the near-independence of 
the measurement fidelity from the relation between the 
charge basis and energy eigenbasis. In Sec. V we discuss 
how decoherence can be introduced to the problem and 
how one should interpret the measurement data in the 
presence of decoherence. We analyze one simple exam- 
ple that shows how the measurement data recorded at 
later times is less 'valuable' than that recorded at earlier 
times. Concluding remarks are given in Sec. VI. Some 
details of the calculations are given in the appendices. 



II. MODEL 

We consider a system composed of a charge qubit ca- 
pacitively coupled to a QPC, as illustrated in Fig. 1. The 
qubit can be viewed as a system where a charged parti- 
cle is trapped in a double-well potential and can occupy, 
and tunnel between, the localized ground states of the 
two wells. We shall denote these states by \L) and 

During the measurement a voltage is applied to the 
QPC, and a current flows through it. Except for the dis- 
cussion given in Sec. V, we shall assume that the QPC 
does not induce any qubit decoherence except that as- 
sociated with the measurement-induced projection. Be- 
fore analyzing the general case, we first consider the case 
where the qubit Hamiltonian is diagonal in the charge 
basis. In this case, there is no mechanism by which the 
states \L) and \R) mix during the system dynamics. As 
a result, if the qubit is initially in the state |L), the long- 
time-averaged QPC current will be given by /l, and the 
qubit will remain in the state \L). A similar statement 
applies to the state \R) of the qubit, with correspond- 
ing QPC current Ir. The QPC current therefore serves 
as an indicator of the qubit 's state in the charge basis 
as long as the qubit Hamiltonian does not 
mix the states of this basis. 

On any finite timescale, there will be fluctuations in 
the QPC current, and the observed value might deviate 
from or I r. The longer the period over which the 
averaging is made, the smaller the fluctuations. One can 
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FIG. 1: (color online) Schematic diagrams of (a) a charge 
qubit with charge states \L) and |-R) measured by a quantum 
point contact (QPC), (b) the probability distributions of the 
possible QPC current values for the two charge states of the 
qubit, (c) the filled energy levels of the QPC's electrodes and 
the energy barriers for tunneling between the electrodes for 
the two charge states of the qubit and (d) a typical QPC cur- 
rent signal. The finite widths of the probability distributions 
in (b) are a result of the finite measurement time. When the 
distance between the two center points, \I r — I l\, is much 
smaller than the widths of the distributions, the QPC per- 
forms a weak measurement on the qubit in the short interval 
under consideration (i.e., the time interval for which these 
probability distributions are relevant). The fast fiuctuations 
in (d) convey the assumption that these fluctuations have a 
much shorter timescale than the qubit's precession period. In 
plotting (b-d) we have assumed that Ir > Ij^, which would 
be the case if the qubit is defined by an extra positive charge 
(e.g. a hole) tunneling between the two wells. Taking the 
opposite case, i.e. Il > Ir, would only result in a few sign 
changes in the analysis below. The main results are not af- 
fected by this choice. 

therefore deflne a measurement timescale that determines 
how long one needs to wait in order to distinguish be- 
tween the states \L) and \R). The relation between this 
timescale and the qubit's Hamiltonian- induced preces- 
sion period separates two measurement regimes: strong 
and weak qubit-detector coupling. As mentioned above, 
this separation is irrelevant when the qubit Hamiltonian 
is diagonal in the charge basis. 

For the remainder of this paper, we analyze the gen- 
eral case where the qubit Hamiltonian is not necessarily 
diagonal in the charge basis. The Hamiltonian of the 
combined qubit 4-QPC system can be expressed as 

H = Hq + Hqpc + -ffint- (1) 

In describing the qubit, we shall use the eigenbasis of the 
qubit Hamiltonian Hq. Thus Hq can be expressed as 

^.= (7|)=-f^. (2) 

where E is the energy splitting between the qubit's two 
energy levels, and a-^ is the 2-axis Pauli matrix. We shall 
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express the ground and excited states of the Hamiltonian 
as |0) and |1), respectively. The electric dipole moment 
operator \R) {R\ — \L) {L\ can be expressed as where 
n is a vector that represents the direction of the charge 
basis relative to the energy eigenbasis: 



<7„ - \R){R\-\L){L\ 
= cos (3 ffz + sin f3 , 



(3) 



where /? represents the angle between the charge basis 
and the energy eigenbasis. The states of the charge basis 
(i.e., the eigenstates of a^) can be expressed as 



\R) = cos^|0>+sin^|l> 
\L) = sin^|0)-cos^|l). 

The QPC Hamiltonian Hqpc can be expressed as 



(4) 



E 

{m, n} G 
tunneling 
channels 



Jn 



(5) 



where the operators am and &„, respectively, annihilate 
electrons in the source and drain electrodes of the QPC, 
alj^ and foj^ are their hermitian conjugates, es,m and eD,n 
are the energies of single-particle quantum states in the 
source and drain electrodes, and Jm,n represents tunnel- 
ing matrix elements between the two electrodes. We take 
the qubit-QPC interaction Hamiltonian to be of the form 



i 



E 

{m, n} £ 
tunneling 
channels 



SJrj 



(6) 



such that the tunneling matrix elements (and therefore 
the QPC current) depend on the state of the qubit in the 
charge basis. The QPC and its operation are illustrated 
in Figs, l(b-d). When J and SJ are independent of m 
and n, the microscopic degrees of freedom of the QPC 
are not affected by, and therefore do not carry any infor- 
mation about, the state of the qubit. This information is 
carried only by the total number of electrons that have 
tunneled through the QPC. In the remainder of this pa- 
per, we shall assume that this is the case. We close this 
section by mentioning that the above description of the 
microscopic structure of the QPC is not needed for the 
analysis below, provided one accepts the description of 
the QPC operation explained in Sec. HI as a phenomeno- 
logical model of the detector. Indeed, the analysis below 
is not restricted to the case of a charge qubit probed by 
a QPC. 



III. MEASUREMENT- AND 
HAMILTONIAN-INDUCED DYNAMICS 

In this section we introduce the theoretical framework 
for analyzing the qubit state evolution and interpreting 
the observed measurement data. 



A. Short-time evolution 

We start our analysis by considering a short time in- 
terval between times t and t + St. For dcfiniteness, we 
assume that during this time interval a large number of 
electrons tunnel through the QPC, such that it is nat- 
ural to define a QPC current I{t) during this short in- 
terval (The current would be defined as the amount of 
charge that has passed through the QPC divided by 5t\ 
The small-current case will be discussed in Sec. IV. F). 
We also assume that a weak-coupling regime exists for a 
properly chosen value of 5t, which means that the QPC- 
current probability distributions (for the states \L) and 
\R)) are broad and almost completely overlap, as shown 
in Fig. 1(b). The coherence time of the QPC and the 
time resolution of the QPC's output signal are both as- 
sumed to be much shorter than the precession period 
of the qubit. As a result, for a properly chosen value 
of 5t one can treat the different time intervals as inde- 
pendent measurement processes, with each measurement 
being made in the charge basis. 

As explained in Appendix A, the QPC-current proba- 
bility distribution for the qubit state \L) can be expressed 
as 



PL[I,5I,8t] = 



2St {siy 

^ ttt™ (A/)' 



■ exp 



2St [I-IlY 



(7) 

and a similar expression can be derived for the qubit state 
\R). Here we have defined a number of useful quantities: 
SI is the size of a finite interval of QPC currents that 
we identify with a single value. One could say that with 
this definition we are turning the probability distribu- 
tions in Fig. 1(b) into histograms with discrete possible 
values for the current I. This definition will be useful for 
both the analytic and numerical calculations below. The 
characteristic measurement time t„i defines the timescale 
needed for the QPC output to distinguish between the 
states \L) and \R), and AI = Ir — II, where for defi- 
niteness we take In > 1^. 

We can now construct matrices that describe the qubit- 
state evolution depending on the observed QPC current 
I{t). When a given value of I{t) is observed in the QPC, 
the density matrix Pq(t) describing the quantum state of 
the qubit is projected (possibly partially) according to 
the observed value. This projection of the qubit's state 
is described by a 2 x 2 matrix (or propagator) that we 
shaU caU UM[Iit),5I,5t]: 
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pS + 6t) = : yUM[I{t),SI,St]p^it)ul[I{t),SI,6t]. (8) 

Tr {[/],,[/(<), <5/, StpMiHt), SI, 5t]pq(t)| 

Measurement theory says that IJM[I{t), SI , 6t] obeys the equation [l7| 

Ul,[lit),51, 5t]UM[I{t), SI, St] = PL[I{t),5I, St] \L) (L\ + P;? [/(t), SI, St] \R) (R\ . (9) 

Since we are keeping only the fundamentaUy necessary terms in the projection dynamics, we ignore any unitary- 
transformation component in UM[I{t), SI , St] that does not affect Eq. ([9]). Any such component can in principle be 
determined experimentally, and it can then be incorporated into the analysis straightforwardly; Such a component 
could result if, for example, the different current states of the QPC have different charge distributions that act 
as an effective gate voltage on the charge qubit. Assuming a large overlap between the QPC-current probability 
distributions, we find that the matrix UM[I{t), SI, St] can be expressed in the form 



UM[Iit),SI,St] = ^PL[I{t),5I,5t] \L) {L\ + ^Pj,[I{t),SI,St] \R) {R\ 



1 



2St (SI) 

TT,r, (AT)\ 



1/4 



/ St[l{t)-I + AI/2]' \ , / St[m~I~AI/2^ 
exp <^ —-2 \ \L) {L\ + exp <^ — \ \R) {R\ 

[ A/ J [ Trn{AT) J 



/ 1,,,,,, St [iit)-7\ 

[ TmAI J [ TmAI 




^AI 



i\R) {R\ - \L) {L\) 



St \i(t) - 71 

= 1+ ^ ^n, (10) 

TmAl 

where we have defined / = (/l + Ir)/^. 

In addition to the measurement-induced evolution described by the matrix Um of Eq. (|f Op . the qubit Hamiltonian 
induces a unitary evolution of the qubit's state. This contribution to the qubit-state evolution (over the time interval 
from t to t + St) is described by the matrix 

UnlSt] = exp i^-iH^Stj « I + i—Sta^. (11) 

The matrices [/^[^(i), SI, St] and UniSt] can now be combined to give the total evolution matrix 

U[Iit),SI,St] w UM[Iit),SI,St]xUH[St] 
St\l{t)-7] 

cx 1+ ^ ^ ^ &„+i-Sta,. (12) 

TmAI 2. 

Note that the high-resolution requirement in the QPC-current signal, which is one of our assumptions mentioned earlier 
in this section, enters in obtaining the above equation: both UM[I,SI,St] and UH[St] are approximately proportional 
to the unit matrix, with lowest-order corrections of order St. Therefore they commute to first order in St, and they 
can be treated as commuting operators when ESt <C 1, i.e. when they both induce small changes to the quantum 
state of the qubit. 

When a given QPC output signal I{t) [starting from the initial time t — Q until the final time tf] is observed, 
one can take the corresponding short-time evolution matrices explained above and use them to construct the total 
evolution matrix ?7Totai[^(^ : — > tf),SI,St]. This matrix describes the evolution of the qubit state from t = until 
t = tf, given that the current signal I{t) was observed (with discretization parameters SI and St). Using the unit 
matrix as the total evolution matrix for t = 0, we find that 

J7Totai[/(i : ^ tf),SI,St] = U[I{tf - St), SI, St] x •• • x U[I{Q),SI ,St]. (13) 



In Sec. III.B we shall explain how measurement information can be extracted from the matrix C/xotai- 
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B. Extracting information from the matrix UTota.i[I{t : tf),5I,St] 



The evolution of the qubit state is obtained by taking into account a long sequence of weak measurements and 
Hamiltonian-induced unitary transformations. This long sequence of events, however, is equivalent to a simple scenario 
in which a single measurement is made on the qubit and the state of the qubit is rotated after the measurement. In 
order to obtain a quantitative description of this alternative interpretation of the measurement data, one takes the 
2x2 matrix C/Totai[-^(^ • — > t/), SI, St] and divides it through to a polar decomposition into two parts, a measurement 
matrix UMcas[I{t ■ — > t/), SI, St] followed by a unitary transformation Uiiot[Iit '■ — *■ i/), SI, St]: 



UTot.i[I{t : ^ tf), SI, St] - URot[I{t : ^ tf),SI, St] x UMo.s[I{t : ^ tf),SI, St] 



(14) 



r 



The matrix UMciis[I{t : ^ t/), SI, St] has the form 

Uucsint : -> tf),SI, St] = y/Fi ItAi) (V'i| + /f^|V^2) 

(15) 

where ]ipi) and ]ip2) are two orthogonal states (the indices 
are assigned such that Pi > P2)- In order to character- 
ize the measurement matrix in terms of the information 
it provides about the state of the qubit, it is useful to 
consider the case where one initially has no information 
about this state. The qubit therefore starts in the maxi- 
mally mixed state 



Pmax. mixed ( q 







(16) 



After a measurement that produces the outcome corre- 
sponding to the matrix UMeas[I(t : tf),SI,St] (and 
ignoring the presence of the rotation matrix UB.ot[I{t : 
tf),SI,St] for purposes of this argument), the qubit 
ends up in a state described by the density matrix 



Paftcr Mcas. 



Pi 



P1+P2 



1^1) ('^iK 



P7 



P1 + P2 



l'/'2) {11^2 



(17) 

This density matrix describes a statistical mixture of the 
states ]'ipi) and IV'2), with a higher probability for the 
former (except for the rare cases where Pi = P2). One 
can therefore say that, as a result of the measurement, 
one now knows that the qubit is more likely to be in 
the state |'0i) than in the state IV'2), in contrast to the 
complete lack of information at the initial time. The 
states ]ipi) and IV'2) therefore represent the measurement 
basis that corresponds to the output signal I{t). The 
measurement fidelity is calculated as follows [ISj] : Let us 
assume that one is given a qubit that is in an unknown 
state, either ]tpi) or ]ip2)- The parameters Pi are the 
probabilities that the outcome defined by I{t), SI and St 
is obtained given that the qubit was initially in the state 
IV'i)- Upon observing this outcome, if one wants to make 
a guess about the state of the qubit, one would maximize 
the probability of making a correct guess by choosing 
the state The difference between the probability 

of actually finding the state and that of finding the 
orthogonal state ]ip2) is given by 



F = 



P^ -P2 



Pi +P2 



(18) 



This expression can therefore naturally be interpreted as 
the measurement fidelity. 

It is worth pausing here to comment on the issue of the 
independence of the measurement basis from the mea- 
surement outcome. One typically thinks of the mea- 
surement basis as being determined by the experimenter 
when designing the experimental setup and the measure- 
ment result (i.e., -1-1 or —1 along the measurement axis) 
being determined stochastically according to the rules 
of quantum mechanics based on the measurement basis 
used in the experiment. The situation considered in this 
paper, however, cannot be described using this simple 
picture. Here we have a two-state quantum system and 
a measurement device that produces one of many pos- 
sible outcomes, not necessarily all providing information 
in the same basis. In this case, the picture of the mea- 
surement basis being independent of the outcome is not 
valid in general. For example, it is possible in a general 
setting for a certain state ]ip) to be a possible result of 
the data analysis presented above, but not the state or- 
thogonal to it. Such a situation arises in the problem 
analyzed in Ref. p^ . and it complicates the interpreta- 
tion of the measurement data. Fortunately, because in 
this paper we are dealing with two symmetric probability 
distributions [see Fig. 1(b)], we find that for every possi- 
ble QPC signal there is an 'opposite' or 'flipped' signal, 
which can be obtained by taking the mirror image about 
the central line in Fig. 1(d) [this central line is defined 
by / = (/l + Ir)/2]. It is not difficult to verify that if a 
given signal corresponds to a certain measured state [?/'), 
the 'opposite' signal will correspond to the orthogonal 
state with the same values of Pi and P2 (see Appendix B 
for the full derivation) . An important consequence of the 
above statement is that the total probability of finding 
one of two opposite signals is independent of the initial 
state of the qubit. Hence, the probability of obtaining 
a certain measurement basis is independent of the ini- 
tial state of the qubit. This fact leads naturally to the 
conceptual interpretation that the measurement basis is 
determined according to some stochastic process, and the 
measurement result along that basis is made after the ba- 
sis has been determined. 

To summarize, the QPC's output signal can be used to 
derive the matrix UTotaiilit : tf), SI , St]. This matrix 
can then be used to determine the measurement basis. 
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the measurement result (i.e., ±1 along the measurement 
axis), the fidehty (or in other words, the degree of cer- 
tainty about the obtained measurement result) and the 
post-measurement rotation. Note that all the analysis 
above is performed independently of the initial state of 
the qubit. If one knows the initial state, one can use the 
evolution matrices in order to determine the probabili- 
ties of different outcomes, as well as the final state of the 
qubit after the measurement. As we shall show below, 
when the measurement fidelity approaches one, the fi- 
nal state can be determined even without any knowledge 
about the initial state. 

It is worth mentioning here that the analysis presented 
above can be applied not only to qubits, but also to sys- 
tems with higher-dimensional Hilbert spaces. The main 
complication in the case of higher dimensions is that the 
description of the results becomes more intricate and less 
transparent. In this context one can contrast the simple 
visualization of any qubit density matrix as a point in 
the so-called Bloch sphere, whereas there is no visually 
simple parametrization for a general density matrix of a 
three-state system. 

C. Incorporating the matrix UTotai[I{t : — » tf),SI,St] 
into the formalism of stochastic master equations 

Equation ^ describes the projection of the qubit's 
density matrix conditioned on the observation of the 



QPC output I{t). That equation treats the evolution 
as occurring in discrete steps. Constructing a master 
equation means taking Eq. ([5]) to first order in 6t, and 
then taking the limit St in order to obtain the corre- 
sponding continuous-time differential equation ^17] . One 
tricky point in carrying out this procedure is the fact that 
the stochastic quantity I{t) — I has fluctuations of size 
1/V6t: 




As a result, when keeping terms only to first order in St, 
one must also keep terms containing {5t)^ x (/(i) — 7) , 
etc. The quantity I{t) — I can be expressed as its ex- 
pectation value at any give point in time, which is given 
by A/(d-„)/2, plus a fluctuation. It is customary to write 
the stochastic master equation not in terms of the current 
fluctuation, but rather in terms of a rescaled quantity, 
SW, whose standard deviation is equal to VSi. It then 
follows that SW should be defined using the relation 

m=7+^{a^) + ^^5W. (20) 

Substituting Eq. into Eq. ^ [ignoring the sub- 
script M here] and keeping terms to first order in St, we 
find that 



Pq{t + St) = Pq - iStHqPq + iStpqHq 



(o'nPq + PqCTn " 2 ((7„) Pq) + (CTnPqCTn " Pq) , (21) 



2 Y^'^m 



where, for compactness, we have not expressed the time dependence of pq explicitly on the right-hand side [In other 
words Pq should be interpreted as Pq{t)]. This equation can now be rewritten in the form 



Pq = 



Hq, Pq 



}'^%- i^nPq + PqO-„ - 2 (o-„) Pq) + {^aPq^n ~ Pq) ■ 



(22) 



r 



where ^(t) is obtained by using the rescaling ^(<) 
SW/St and taking the limit St — 0: 



mm) = 5{t-t'), (23) 

where S{t — t') is the Dirac delta function. Equation (1^^ 
is the stochastic master equation for the system under 
consideration. 

The question now is how the matrix ?7Totai[-^(^ ■ — > 
tf),SI,St] fits into the picture of the stochastic master 
equation. One can derive a stochastic equation for the 
matrix C/Totai[-^(^ '■ — > tf),SI,St] by considering the 
change caused by a matrix of the form given in Eq. ; dtl- 
note that the evolution of the matrix UTotai[I{t : ^ dt 



tf),SI,St] is governed by Eq. p^ . One problem that 
arises here is that the matrix in the first line of Eq. ([72]) 
tends to zero when St — > 0. Taking this limit therefore 
leads to ill-defined quantities. In order to avoid these 
problems, we deal with a re-normalized version of the ma- 
trix U[I{t),SI, St], as given in the second line of Eq. (fT2|l : 
it should be noted here that normalization of the matrix 
C^Totai is irrelevant for purposes of extracting the mea- 
surement basis and fidelity. It is now straightforward to 
see that an unnormalized version of the matrix t/xotai 
obeys the equation 



Total 



rotal- 



(24) 
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This equation supplements the stochastic master equa- 
tion for purposes of interpreting the QPC's output signal 
in terms of a measurement result. In this context, one 
should note that what Eq. gives is the qubit's den- 
sity matrix as a function of time. When /3 7^ and the 
charge states mix during the measurement, the density 
matrix at later times can become completely different 
from that at the start of the measurement process. As- 
suming that the experimenter does not have any informa- 
tion about the qubit's state at t = 0, he would be able to 
use Eq. to determine the qubit's state at later times, 
but he would not be able to cast the results in terms of 
information that he has learned about the qubit's state 
at i = 0. The matrix L/xotai, through the decomposition 
explained in Sec. III.B contains the information needed 
to make such a statement about the measurement result. 
In fact, one can say that, when taken in combination 
with the initial qubit density matrix, Eq. (j24p replaces 
Eq. (|22p : in addition to containing the measurement in- 
formation, the matrix [/xotai can be used to calculate the 
density matrix at any time. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We now present the results of our numerical calcu- 
lations. Typical parameters of the numerical calcula- 
tions are as follows: The discrete steps are taken to be 
St = O.Itt/E, such that one period of the Hamiltonian- 
induced coherent oscillations would be divided into 20 
steps. The QPC current probability distributions are dis- 
cretized into 100 possible current values: each distribu- 
tion is Gaussian with a standard deviation a of 10 steps 
(thus there are about 50 possible values for the QPC 
current with non- negligible occurrence probability). The 
distance between the two center points of the current 
probability distributions A/ determines the qubit- QPC 
coupling strength. This parameter is varied when analyz- 
ing the effect of the qubit-QPC coupling strength on the 
behaviour of the system. The number of time steps in a 
single calculation was different for different calculations. 
However, this number was always large enough that the 
final qubit state reached perfect purity up to less than 
one part in 10^, ensuring that none of our numerical re- 
sults were sensitive to the exact choice of this parameter. 
Each measurement procedure is repeated up to 10^ times 
in order to obtain accurate statistical averages. We have 
also used parameters different from the typical ones given 
above and found that the results are unaffected by the 
specific choice of these parameters. 

A convenient parameter for purposes of characterizing 
the qubit-QPC coupling strength is Et„i/ (27r), where t„i 
is the timescale needed to obtain sufficient QPC signal 
to read out the state of the qubit (for the time being one 
can think of this definition as applying to the case when 
/3 = 0; but see below). As the standard deviation of the 
QPC-output effective probability distribution scales as 
a / \/N when the measurement step is repeated TV times. 



the measurement time Tm can be naturally defined as the 
product of the time step 5t and the value of N at which 
2a/^/N = A[. T„i is therefore given by 



We shall use the parameter £'T,„/(27r) for quantifying 
the qubit-QPC coupling strength, and we shall use as 
a characteristic measurement timescale when presenting 
our results below. 



A. Stochastically determined measurement basis 
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FIG. 2: (color online) The spherical coordinates 6 and cj> 
defining the stochastically determined measurement bases ob- 
tained in simulations of the experiment under consideration 
(note that, since the measurement basis is defined by an axis 
extending in two opposite directions, only the upper hemi- 
sphere is needed; any measurement axis will either he in 
the plane of the equator or have one end above the equa- 
tor). Each figure contains 200 points. In each row, the 
qubit-QPC coupling strength is kept fixed and /3 is varied: 
/3 = 0, 7r/8, 7r/4, Stt/S and 7r/2. In each column /3 is kept 
fixed and the qubit-QPC coupling strength is varied. In the 
top row, ETm/Cin) = 0.01, which is deep in the strong- 
coupling regime; The measurement basis is very close to the 
charge basis in about 99% of the runs. In the middle row, 
ETm/(2Tr) — 0.2, which can be identified as the intermediate- 
coupling regime; The measurement basis deviates substan- 
tially from the charge basis in about one half of the runs. In 
the bottom row Erm/ {2tt) — 5, which is in the weak-coupling 
regime; The measurement bases are, in general, spread over 
the entire hemisphere. In generating this figure, the initial 
state was taken to be the maximally mixed state (note that 
the results are independent of the initial state). 

First, in Fig. 2 we show the spherical coordinates 9 
and of a large number of (stochastically determined) 
measurement bases for different qubit-detector coupling 
strengths and different values of (3. In the strong-coupling 
regime (top row in Fig. 2), the parameter Erm /{in) is 
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much smaller than one, and the measurement is com- 
pleted before any Hamiltonian-induced dynamics can oc- 
cur. As a result, the measurement basis is always the 
charge basis, i.e. the natural measurement basis for the 
detector under consideration [note here that the charge 
basis is characterized by the direction {9, (p) ~ {(3, 0)]. As 
the qubit-detector coupling strength is reduced (middle 
row in Fig. 2), the measurement bases start to deviate 
from the charge basis, and they develop some statisti- 
cal spread. This region could be called the intermediate- 
coupling regime. Deep in the weak-coupling regime (bot- 
tom row in Fig. 2), the measurement bases are spread 
over all the possible directions. The probability distribu- 
tion of possible measurement bases will be analyzed in 
Sec. IV.D. 

When f3 — (leftmost column in Fig. 2), the mea- 
surement basis is always the charge basis, regardless of 
the qubit-QPC coupling strength. This result is natural, 
since this case is the simple one with no mixing between 
the states \L) and \R). When (3 = 7r/2, all the possible 
measurement bases lie in the x-y plane. The reason be- 
hind this result lies in the fact that this situation is equiv- 
alent to one where one makes a large number of weak 
measurements in the x-y plane with no Hamiltonian- 
induced precession. As explained in Appendix C, the 
resulting measurement bases can only be in the same 
plane as the actually performed measurements, which is 
the x-y plane in this case. 

The fact that the measurement basis is generally un- 
predictable, and therefore uncontrollable, is a rather 
strange phenomenon from a fundamental point of view. 
From a practical point of view, one can wonder whether 
anything useful can be done with such measurements that 
are performed in a stochastically determined basis. If one 
absolutely requires a measurement in a given basis, mea- 
surement results in different bases would be less useful. 
One could then treat the deviation of the observed mea- 
surement basis from the desired one as an experimental 
error and deal with it accordingly. This point will be 
discussed further in Sec. IV. E. 



1.0 
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FIG. 3: (color online) The measurement fidelity as a func- 
tion of measurement duration for three different values of the 
angle /3 between the charge basis and the energy eigenbasis: 
l3 — (red; lowest line), 7r/4 (green) and tt/2 (blue). Here 
ETm/Cin) — 5, i.e. deep in the weak-coupling regime. The fi- 
delity increases from zero to one as the measurement duration 
increases, and it is almost independent of the angle /3. 



extract useful measurement information when [3^0, the 
information acquisition rate is not reduced, and in fact 
enhanced, by the Hamiltonian-induced precession. 

In the limit when the fidelity reaches one, the matrix 
UMcas[I{t ■ tf),SI,St] has the form \^pi) (-0i|. It is 
straightforward to verify that any further evolution of the 
system will not aff'ect UMeas[I{t ■ tf),6I,6t]. In this 
case, further measurement only changes the final state 
of the qubit, which would be given by t/Rot[-^(^ '■ — > 
tf),6I,6t]\iji). 

It is worth mentioning here that the fidelity for the 
case f3 — is very well fitted by the function 



B. Information acquisition rate 

In Fig. 3, we plot the measurement fidelity as a func- 
tion of measurement duration for three different values 
of P (i.e., the angle between the charge basis and the 
qubit 's energy eigenbasis), keeping all other parameters 
fixed. We can see that the fidelity approaches one for 
long enough measurement time, regardless of the value 
of the angle (3. Furthermore, the fidelity is almost in- 
dependent of p. In fact, and rather counter-intuitively, 
it turns out that the fidelity for /? is higher than 
that for P — 0. The reason behind this phenomenon is 
that one can obtain a higher fidelity by performing weak 
measurements along different axes rather than along the 
same axis, as explained in Appendix C. This result shows 
that even though more complicated analysis is needed to 



where erf stands for the error function. This result can 
be obtained using the following calculation. We take the 
probability distribution in Eq. (|57D and we replace St by 
t. We now identify the portion of this distribution that 
corresponds to / > {II + Ir)/'2 as being the probabil- 
ity of making a wrong inference about the state of the 
qubit. The difference between the probability of making 
a correct inference about the state of the qubit and that 
of making a wrong inference gives Eq. . The fidelity 
in Fig. 3 is also relatively well fitted by the function 

F(t)«l-exp(- — 1, (27) 

L Trn J 

which is an indication that our definition for the mea- 
surement time Tm (see Appendix A) is a reasonable one. 
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C. The information contained in the QPC's output 

signal 

The QPC produces an output signal that contains fast 
and large fluctuations. There are therefore a very large 
number of possible output signals that the QPC can dis- 
play. We now ask the question of what information is 
contained in a given QPC output signal. It turns out 
that only a small amount of information in the signal 
concerns the state of the qubit. As explained above, one 
can follow a conceptually straightforward calculation in 
order to extract the measurement basis, measurement re- 
sult and fidelity from the noisy signal. The signal also 
contains information about how the qubit 's state is ro- 
tated after the effective measurement. It should be noted 
here that measurement bases and post-measurement ro- 
tations generally come in all possible combinations. The 
measurement basis, result and fidelity on the one hand 
and the post-measurement rotation on the other hand 
can therefore be seen as two separate pieces of informa- 
tion that are contained in the QPC's output signal. The 
rest of the information contained in the highly noisy sig- 
nal, which is the vast majority of information contained 
in the signal, concerns the QPC itself. Thus the large 
fluctuations that one observes provide information about 
the state of the QPC in the specific experimental run un- 
der consideration. For purposes of reading out the state 
of the qubit, this information is discarded. This situ- 
ation is most clearly demonstrated in the case (3 = 0, 
where one is only interested in the time average of the 
signal, discarding all the fluctuations around the average. 





D. Spectral analysis of the measurement-basis 
probability distribution 

We now take sets of stochastically determined mea- 
surement bases (i.e., results of the kind presented in 
Fig. 2) and use them to extract probability distribu- 
tions 5(0, (/>) for the measurement basis. The probability 
distributions arc normalized to unity: J g{9, (j))dD, = 1 
with the integral covering the hemisphere < < 7r/2, 
< < 27r. In particular, we are interested in these 
probability distributions in the weak-coupling regime. In 
this regime, the probability distributions will be symmet- 
ric about the z axis, i.e. g{9, (j)) will be independent of 4>. 
We use this fact and write the probability distribution as 



oc 



71=0 



(28) 



FIG. 4: (color online) Spectral analysis of the probability dis- 
tribution g{6,(f>): (a) The spectral components c„ for n = 
(red circles), 2 (black squares), 4 (green x symbols), 6 (blue 
triangles), 8 (purple + signs) and 10 (green diamonds) as 
functions of _ETm/(27r) for j3 = tt/8. Note the logarithmic 
scale on the x axis. As Erm/ (Stt) is increased, the spec- 
tral components clearly reach very stable values by the time 
ETm/{2n) reaches one. These results indicate that the value 
ETm /{27r) = 5 is sufficiently deep in the weak-coupling regime 
for purposes of studying the limit J5rm/(27r) — > 00. (b) The 
spectral components Cn for n = 0,2,4 and 6 [same symbol 
usage as in (a)] as functions of f3 in the weak-coupling regime 
[here we use the value Erm/{2n) = 5]. As /3 is increased, the 
probability distribution g{6,cf>) gradually changes from being 
concentrated at the energy eigenbasis to being spread over all 
directions and then to being concentrated on the equator per- 
pendicular to the energy eigenbasis. There seems to be one 
value of (3 [close to O.Svr, which is close to cos~^(l/\/3)] where 
all the spectral components with n 7^ vanish, i.e. g{6, (j)) be- 
comes completely uniform. 



where F„_o( 
by 



are the usual spherical harmonics given 



2n+ 1 
47r 



F„(cos( 



(29) 



and Pn{x) are the Legendre Polynomials. Note that since 
we only need to consider the hemisphere < 9 < 7r/2, 



we only need to keep even values of n in the above series. 
The coefficients c„ are given by 



Cn = 2 d9sin9 d(j,g{9,(j))Y^^o{9A)- (30) 
Jo Jo 



10 



In our calculations, we have a collection of N stochasti- 
cally determined measurement bases. The coefficients c„ 
can be calculated from these data sets using the formula 

^n^^Y.^nA(^J,^j)- (31) 

j 

where the index j labels the different data points. 

In order to ensure that we are using parameters that 
are sufficiently deep in the weak-coupling regime, we start 
by performing the above spectrum analysis while varying 
i?Tm/(27r). The results of this calculation for /3 = vr/S 
are shown in Fig. 4(a). The coefficients c„ seem to be 
almost independent of i?T„i /(27r) when this parameter 
is larger than one. From these results we conclude that 
the point i?T„j/(27r) = 5 is sufficiently deep in the weak- 
coupling regime for purposes of calculating g{0, (p) in the 
limit i?T,„/(27r) — > oo. 

We now calculate c„ (up to n — 10) in the weak- 
coupling regime for different values of (3 and plot the 
results in Fig. 4(b). The probability distribution changes 



gradually as we go from /3 = 0, where the measurement 
basis is always along the z axis (here the charge basis) , to 
(3 = it/2, where the measurement basis is always in the 
x-y plane. There is one point (/? « O.Stt) where g{9, (f>) is 
either approximately or exactly uniform over the entire 
hemisphere. 



E. Short-time dynamics 

We now consider the case where the measurement du- 
ration is short enough that the measurement fidelity is 
much smaller than one. In this case the state of the qubit 
experiences only a small amount of projective evolution. 
As a result we can calculate analytically the total effect 
of the projection and coherent dynamics. 

Dividing the total time interval under consideration 
into a large number of infinitesimal intervals, we can 
write the total evolution matrix as 



f/Totai [I{t : -> t/), SI, 5t] « [1 + <5p(i/)<7„] e-^^**/2 . . . e-^^^*/^ [l + Sp{St)&n] e'*^^^*/^ ^ Sp{Q)an] 

= e^-^s*//2 [1 + <5p(t/)a„(t^)] X ... X [1 + dpiSt)an(st)] x [l + 5p(0)a„(o)] 
« e^'^'^^f /'^ [1 + Spit f) {sin 13 coa{Et f)a.^ + sin (3 sin{Et f)ay + cos (3&^}] x ... 

X [1 + Sp{0) {sin/3cos(0)(J:r + sin/3sin(0)(Tj^ + cosP&z}] 

« e**--Et//2 f ' dt^^ {sin (3 cos{Et)aa,+ sin (3 sin{Et)ay + cos (3&z} 

Jo St 

= e"'^''''/^[l+P^<J^+Pyay+p,a,], (32) 



where we have defined 

,,,, = ^IM^ 

n{t) = sin/3cos(£'t)a; -f sin/3sin(_Ei)y -I- cos/3z 
Px = sin /3 / — cos{Et)dt 



Py = sin f3 



T,„A/ 



I a T,nAI 

Pz — cos/3 / —dt. 

Jo TmA/ 



^ sin{Et)dt 



(33) 



Taking the total measurement duration tf to be a mul- 
tiple of the qubit's free-precession period, the different 
measurements strengths Px,y,z are different Fourier com- 
ponents of the observed QPC signal. These Fourier com- 
ponents can be thought of as independent quantities, 
leading one to think of the net result as three independent 
weak measurements, one along each of the three axes x, 
y and z. 



We have performed spectral analysis as was done in 
Sec. IV. D on the possible measurement outcomes in the 
short-time case analyzed above (we do not show the full 
results here). In the calculations px, Py and pz were as- 
sumed to be Gaussian-distributed quantities with root- 
mean-square values that obey the relation 



rms . „rms . rms 
-Pa; . -Py . -Ps 



sin (3 sin /3 



cos (3. 



(34) 



The overall features of the dependence of c„ on f3 are 
similar to those seen in Fig. 4(b), but there is no exact 
agreement between the two cases. For example, the coef- 
ficient C2 is close to being a straight line all the way from 
P ~ to f3 — 7r/2 in the present case, as opposed to the 
curved path that the black squares follow in Fig. 4(b). 
Finally, we ask whether there is a value of /3 at which all 
three measurements have the same strength on average. 
Using Eq. it is straightforward to see that this situa- 
tion occurs when cos (3—1 /VS, i.e. /3 = 0.3047r; the point 
at which all measurement directions are equally probable 
in Fig. 4(b). The calculations of this subsection therefore 
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give us a simple interpretation of the significance of the 
point (3 = O.Stt in the present context. 



F. Alternative interpretation: Choosing the 
measurement basis independently of the 
measurement outcome 



The idea that the measurement basis is determined 
stochastically and is only specified when the experiment 
is completed is conceptually unsettling. We therefore 
take a different approach in this section. We now re- 
quire that in order to have a meaningful measurement, 
the measiircmcnt basis must be specified by the setup 
and not by the stochastically determined measurement 
outcome. 

In the strong-coupling regime, there are no conceptual 
problems. In this case, the measurement is almost in- 
stantaneous, and the measurement basis is clearly the 
charge basis. The measurement fidelity is almost 100%, 
up to errors caused by the finiteness of the measurement 
time compared to the qubit-precession period. These er- 
rors can be quantified straightforwardly, as we shall do 
shortly. 

The weak-coupling regime is the one where the concep- 
tual problems arise. If we take the case where the distri- 
bution of measurement bases is essentially uniform over 
all direction (which seems to be the case when /3 w O.Stt), 
we again obtain a simple, yet somewhat surprising, re- 
sult. By symmetry we can argue that all directions are 
equivalent. We can therefore choose any basis as the 
measurement basis. Based on this choice, we then di- 
vide all the different possible outcomes (i.e. the different 
possible QPC output signals) into two sets, one of which 
corresponds to the measurement result +1 and the other 
of which corresponds to the result —1 in the chosen ba- 
sis. An important question that arises in this scenario is 
what the average measurement fidelity would be under 
such 'worst-case' conditions of having no preferred mea- 
surement basis. In the following, we shall denote this 
average measurement fidelity by the symbol F in order 
to stress the difference between the results of this subsec- 
tion and those of Sec. IV. B. For the completely uniform 
case, the average fidelity F can be evaluated as follows: 
let us take a spin pointing in the +z direction and assume 
that we have specified the z axis as the one defining the 
measurement basis. If we measure the spin along an axis 
that makes an angle 6 with the z axis, we obtain the re- 
sult -hi with probability cos^(0/2) and the result —1 with 
probability siv?{9/2). By dividing the possible outcomes 
into two sets (with all outcomes that correspond to states 
in the upper hemisphere being interpreted as -|-1 along 
the z axis and all outcomes that correspond to states in 
the lower hemisphere being interpreted as — 1 along the z 
axis) and taking the difference between the probability of 
correctly obtaining the result +1 and the probability of 
obtaining the wrong result —1, we find that the average 



fidelity (averaged uniformly over the entire hemisphere) 



uniform, Er^ri — >oo 



tt/2 



I 2 ^ 9 ^ 

dO sin 9 cos - — sin - 



(35) 



The fidelity in this 'worst-case scenario' is therefore 50%. 

We emphasize again that this fidelity is independent of 
our choice of measurement basis. This result again con- 
trasts the usual intuition where the outcomes of mea- 
surements performed in one basis provide no information 
at all in an orthogonal basis. In the present case, all 
measurement bases are equivalent, and one would obtain 
the same (average) measurement fidelity in any one of 
them. This result fits well with the result that will be 
explained in Sec. IV. G that quantum state tomography 
can be performed by repeating the same measurement 
procedure a large number of times. This measurement 
procedure gives the same amount of information in all 
different bases. 

As explained above, the lowest possible average fidelity 
is 50%, and it occurs when (3 « O.Stt. Two other values of 
/? allow for simple reasoning. The obvious one is the case 
of /3 = 0, where the measiucmcnt basis is {\L) , \R)} and 
the fidelity is 100% for all coupling strengths. The other 
exception occurs when /3 = tt/2. As can be seen from 
the rightmost column of Fig. 2, the measurement basis is 
always in the x-y plane. In the weak-coupling regime, the 
possible measurement bases are spread uniformly around 
the equator. In that limit, one can choose any axis in the 
x-y plane as defining the measurement basis. All these 
bases give the same value for the average fidelity: 



1 /-Stt 



2^ 
COS - 



sin^- 



(36) 



Numerical results for the average fidelity for differ- 
ent values of (3 in the weak-coupling limit are plotted 
in Fig. 5. First, in Fig. 5(a) we plot the average fidelity 
as a function of the chosen measurement angle ^Meas for 
different values of /3. We find that for f3 < O.Stt choosing 
the energy eigenbasis as the measurement basis gives the 
highest average fidelity, whereas for (3 > O.Stt choosing 
a basis that is orthogonal to the energy eigenbasis gives 
the highest average fidelity (note that there is an infinite 
number of such bases). In Fig. 5(b) we plot the average 
fidelity as a function of /3 for three choices of the mea- 
surement basis: ^Meas = 0, ^Meas = tt/2 and the charge 
basis. It should be noted here that in the weak-coupling 
limit there is symmetry about the z axis. As a result, 
any axis that is obtained by rotating the charge basis 
about the Hamiltonian axis will result in the same value 
of the fidelity as the charge basis. Choosing the charge 
basis as the measurement basis never results in an aver- 
age fidelity that is substantially lower than the maximum 
possible value. 
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In the intermediate coupling regime, one can choose 
the measurement basis by maximizing the average fi- 
dehty. For example, for the parameters of the middle 
row in Fig. 2 one can see that the measurement basis 
that would maximize the fidelity deviates slightly from 
the charge basis (unless (3 — 0). One can understand 
this result naturally by keeping in mind that the qubit's 
state precesses about the Hamiltonian's axis while being 
measured. The measurement basis is therefore approxi- 
mately obtained by taking the charge basis and rotating 
it by the appropriate angle about the qubit Hamiltonian 
(in the opposite sense from state precession). 

One could perform a numerical calculation in order 
to analyze the behaviour of the maximum-fidelity mea- 
surement basis as the qubit-detector coupling strength 
is varied. However, here we focus on a question that 
seems more experimentally relevant, namely analyzing 
the measurement fidelity with the charge basis chosen as 
the measurement basis. This fidelity is shown in Fig. 6 
for the two cases /3 — 7r/4 and /3 = tt/2. In Fig. 7, 
we plot the absolute value of the quantity 2TTdF/d{ETm) 
[i.e., as in the initial slope in Fig. 6] as a function of 
/3. This quantity can be used to obtain an estimate for 
the measurement errors associated with having a finite 
measurement time. 



G. Small-current case 
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In the above analysis, we have assumed that the QPC's 
output is essentially a continuous variable with a large 
number of possible values for the current at a given in- 
stant. If the QPC's current is so small that on the 
timescale of monitoring the output signal very few elec- 
trons can tunnel through the QPC, the physical picture 
changes substantially (here we assume that the number 
of electrons that have tunneled through the QPC can 
be measured accurately). One now has a small num- 
ber of possible values for the QPC's output signal. The 
discretization used in our analysis, which was done for 
computational convenience above, is now provided natu- 
rally by the design of the setup. One can then follow a 
similar analysis to that of Sec. Ill concerning the grad- 
ual progression of the measurement and the evolution of 
the quantum state of the qubit. Note that considering 
the small-current case only makes sense if the number 
of electrons that tunnel through the QPC during a full 
qubit-precession period is small. Otherwise, one can al- 
ways choose between the small- and large-current cases 
through the choice of 8t. Clearly, the exact choice of 
parameters used in the theoretical analysis cannot affect 
the physical results. 

An important difference between the cases of small and 
large currents is that one expects the probability distribu- 
tions shown in Fig. 1(b) to have almost symmetric shapes 
for the large-current case, whereas they almost certainly 
are not symmetric for the small-current case. This point 
is most clearly seen by considering the situation where 



FIG. 5: (color online) The average fidelity F when defining 
the measurement basis independently of the specific (stochas- 
tic) outcome in the weak-coupling regime [i5rm /(27r) = 5]: 
(a) The average fidelity as a function of the 'manually' cho- 
sen measurement angle 0Mcas for five difi'erent values of /3: 
(red circles), vr/S (blue squares), 7r/4 (green triangles), Stt/S 
(black -I- signs) and 7r/2 (cyan x symbols). Note that the 
results are independent of the choice of the azimuthal an- 
gle </>Meas, which IS also chosen manually along with ^Moas. 
One can see that for any value of /3 the maximum fidelity 
is obtained by choosing either ^Moas = (for /3 < O.Stt) or 
^iMeas = 7r/2 (for /3 > O.Stt). For (i = 0.3 (not shown) all 
values of ^Mcas give essentially the same value for the fidelity, 
namely 50%. (b) The average fidelity as a function of /3 for 
three different choices of measurement angle ^Mcas: ^Moas = 
(blue squares), 0Moas = l3 (red circles) and ^Moas = n/2 (green 
triangles). The choice ^Moas ~ P always gives a near or true 
maximum in the fidelity. However, with the exception of a 
single point around j3 = O.Stt, the maximum occurs at either 

^Mcas = or ^Moas = TT /2. 



at most one electron tunnels though the QPC during the 
short interval under consideration. The analysis and re- 
sults in this case follow closely the analysis and results for 
a switching-based detector, studied recently in Ref. p^ . 
The fact that, unlike the case of a switching-based detec- 
tor, the measurement now continues after the occurrence 
of a tunneling event can be incorporated into the analy- 
sis straightforwardly. Our argument concluding that the 
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FIG. 6: (color online) The average fidelity F as a function 
of ETm /{2Tr) when the measurement basis is decided (by the 
experimenter) to be the charge basis. The blue open circles 
correspond to the case /3 — 7r/4, and the red filled circles cor- 
respond to the case /3 — Tr/2. The dashed lines represent the 
asymptotic values of the fidelity in the weak-coupling limit: 
0.54 for the case /3 — 7r/4 and 2/tt for the case /3 = n/2. 




FIG. 7: The rate of charge-basis fidelity deterioration with 
increasing measurement time \2TTdF / d{ETm)\ as a function of 
l3 in the strong-coupling regime. In other words, the quantity 
on the y axis is the absolute value of the slope of curves similar 
to those shown in Fig. 6 at Erm/2'K = 0. Each data point is 
determined by calculating the fidelity F at ETm/{2n) = 0.06 
for a given value of /3 and then taking the ratio [1 — F{x)]/x, 
where x = ETm./{2-K). 



choice of the measurement basis is independent of the 
initial state of the qubit breaks down. Perhaps more im- 
portantly, the conceptual picture that the measurement 
basis is chosen first (according to some probability dis- 
tribution) and the measurement result is then obtained 
in that basis becomes invalid. For example, it is now 
possible for a certain state to be a possible result 
of the data anal ysis presented above, but not the state 
orthogonal to it p^ . 



It is also worth mentioning in this context the case 
where the tunneling process is coherent on the timescale 
of qubit precession. In this case the QPC measures the 
qubit in the energy eigenbasis, regardless of the angle /3 
(except when P = tt/2 where the QPC fails to perform 
any measurement on the qubit). For further discussion 
of this case, see e.g. [5i], and for a discussion of possi- 
ble advantages of weak measurement when dealing with 
multi-qubit circuits, see e.g. [loj . 



H. Quantum state tomography 

One example of a procedure where the uncontroUabil- 
ity of the measurement basis can be harmless is quan- 
tum state tomography (QST). In fact, the original idea 
of QST included performing measurements in all pos- 
sible bases [1^. We have simulated QST by repeating 
the measurement procedure a large number of times, 
obtaining a set of measurement results (in the form of 
pre- measurement qubit states) , and then minimizing the 
function (see Appendix D) 



T(r, 0, ^) = ^ [1 - r cosl](0, 0, 0„ </.,)]' 



(37) 



where r, 6 and 4> are the spherical coordinates of a point 
in the Bloch sphere; j is an index labelling the differ- 
ent runs of the experiment, the direction {9j,(j)j) defines 
the qubit state obtained in a given measurement, and 
il{d,(l),9j,(j)j) is the angle between the directions {0,(j)) 



and 



). We have chosen several initial states cover- 



ing the Bloch sphere, and the tomography procedure al- 
ways produced correctly the initial state of the qubit. For 
the largest value of the qubit-detector coupling strength 
that we used [see top row of Fig. 2], the procedure be- 
came unreliable, because the vast majority of the mea- 
surements are performed in one basis. 

It is worth pausing here to ask the question of the min- 
imum requirements for a single measurement setting to 
provide full QST information. The qubit's density ma- 
trix is characterized by three parameters. One therefore 
needs the ensemble of measurements to produce three in- 
dependent variables. If the detector has four possible out- 
comes (or output signals), one would obtain three inde- 
pendent probabilities for three of the outcomes (with the 
probability of obtaining the fourth outcome determined 
by probability normalization). A four-outcome detector 
is therefore sufficient, in principle, to perform full QST 
on a qubit. One can generalize this argument straight- 
forwardly to the case of larger systems: One needs a 
detector with at least as many possible outcomes as the 
square of the size of the Hilbert space in order to perform 
QST in that Hilbert space (we do not necessarily imply 
pure states here). 

In Ref. an 'all-in-one' measurement procedure was 
analyzed, where the measurement results of identically 
prepared setups (using a switching-based detector) can 
be used to perform full state tomography of the initial 
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qubit state as well as extract all relevant system param- 
eters. One could ask whether a similar situation occurs 
in the qubit-QPC system. The answer is yes, in prin- 
ciple. However, because the analysis of Ref. was 
almost entirely analytic, it was possible to obtain sim- 
ple expressions to which the measurement data would 
be fitted in order to extract all the available parameters. 
In this paper, we have not been able to obtain similar 
analytic expressions. As a result, it would require exten- 
sive numerical calculations to compare the observed data 
with those that would be expected for all possible sets of 
parameters and then find the parameters that produce 
the best fit. Given the relatively long computation times 
required for even the simple numerical calculations pre- 
sented in this paper, we have not attempted to simulate 
the full parameter-extraction procedure. 

V. DECOHERENCE 

So far we have assumed that there is no decoher- 
ence in the system other than that associated with the 
measurement-induced projection. In this section we 



briefly consider the effect of additional decoherence chan- 
nels on the measurement process. 

As a first step, one can make the following observa- 
tions: the measurement has a characteristic time Tm, 
which determines how long it takes for the experimenter 
to obtain a high-fidelity measurement result. Decoher- 
ence introduces its own decoherence timescales, which 
we collectively denote by r^. If Tm ^ t^, decoherence 
should have little effect on the analysis and results of the 
previous sections. If, on the other hand, ^ t^j, de- 
coherence will distort the state of the qubit before any 
reliable measurement result is obtained. This situation 
renders the measurement almost completely useless. One 
should keep in mind, however, that the above statement 
depends strongly on the nature of the decoherence chan- 
nel. For example, if the energy eigenbasis coincides with 
the charge basis (i.e., there is no coherent dynamics mix- 
ing the states \L) and \R)) and there is no relaxation 
between the states \L) and \R), pure dephasing between 
these states will not hamper the measurement, regardless 
of the dephasing time. 

Decoherence can also be introduced to the stochastic 
master equation such that Eq. (|22p turns into: 



H, 



q; Pq 



m 



((TnPq + PqtTn 



2 (o-n) Pq.) + -. (CTnPqCTn 



Pq)+Z^7/c 



m 



where the index k represents the different decoherence 
channels, with 7^ and Ak being the rate and quantum 
operator that describe a given decoherence channel. 

Modifying Eq. ([24]) in order to include decoherence is 
somewhat trickier. The difficulty lies in the fact that the 
evolution of the density matrix is no longer of the form 

Pq ^ UpqU\ (39) 

ignoring here terms in Eq. (j22p that were added for nor- 
malization purposes. One way to deal with this new dif- 
ficulty is to note that the evolution of the density matrix 
caused by decoherence is linear and therefore write the 
qubit's density matrix in vector form, i.e. the transpose 
of {poo, Pii: Poi: Pio)- With this rearrangement of the ma- 
trix elements, the evolution matrices are now expressed 
as 4 X 4 matrices acting on the density matrix (which now 
looks like a vector) from the left: 

p^~*Upq, (40) 

where we now use the symbol U in order to distinguish 
the 4x4 matrices of this section from the 2x2 matri- 
ces that are applied symmetrically on both sides of the 
density matrix. 

We shall not write down the detailed expressions for 
the evolution matrices U, because they would not add 



insight to the present discussion. Instead we point out 
that deriving a master equation for these matrices can 
be done straightforwardly by rearranging Eq. (|38p into 
vector format. The differential equation that describes 
the evolution of the density matrix pq also describes the 
evolution of the propagator Z^Totai , similarly to the fact 
that the Schrodinger equation can be applied to wave 
functions or propagators. We now assume that the total 
evolution matrix Z^xotai (including the effects of measure- 
ment, Hamiltonian-induced dynamics and decoherence) 
has been obtained using such a stochastic equation, and 
we ask the question of how to extract measurement in- 
formation from it. The answer is a generalized version 
of the one obtained in the absence of decoherence. The 
total evolution matrix can be decomposed into two matri- 
ces: a measurement matrix followed by a trace-preserving 
'quantum operation'. This quantum-operation part de- 
scribes both unitary evolution and decoherence [231 ■ The 
measurement matrix can be used to extract the measure- 
ment basis, result and fidelity, just as was done in the 
absence of decoherence. 

We now give a simple example illustrating the inter- 
play between measurement and decoherence. We assume 
that the Hamiltonian is diagonal in the charge basis, 
and therefore there is no coherent dynamics between the 
states \L) and \R). A measurement matrix that corre- 
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spends to an outcome that occurs with probability P for 
the maximally mixed state and has fidelity M now takes 
the form 



( l±M 



V 














V'l-M2 




\ 



VI - M2 / 

(41) 

A relaxation matrix with relaxation factor Dj. (which can 
be thought of as decaying from one to zero as time goes 
by) is given by 



^relaxation (-^^r) — 



/1 1 


- Dr 








\ 





Dr 























\o 








^/D'r 


J 



(42) 



and a dephasing matrix with dephasing factor Dp is given 

by 



^dephasing (-^p) 




(43) 



If a qubit is measured, with the measurement outcome 
corresponding to the matrix UueasiP, M), and the qubit 
then relaxes and dephases according to the matrices 
Wrciaxation(-Dr) and Z^dephasing(-Dj,), the total evolution 
matrix describing this sequence of events is given by 



^Total 



= p 



^dephasing (-Dp)ZYrelaxation {Dr )^Meas (-P) -^) 

/ 1±M {l-Dr){lTM) 
Dr{lTM) 





V 



M2 







Dr, 







/D^y/T 



(44) 



M2 J 



This matrix will be used as a reference matrix for the scenario that we describe below. In particular, we shall use the 
fact that the measurement fidelity is given in terms of the matrix elements in the top left corner: 



F = 



U 



1,1 

Total 



1,2 
Total 



-Ur 



2,2 
Total 



u 



Total 



■ u 



Total 



-u 



Total 



(45) 



We now consider a qubit that is constantly being probed by the QPC while undergoing gradual relaxation and 
dephasing. The measurement-induced evolution is now expressed as 
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(46) 



(47) 



Using the above form for the evolution matrix and similar ones for relaxation and dephasing, we arrive at the equation 
of motion for (an unnormalized version of) the total evolution matrix: 



^^Total = St 



(G{t) 



V 



-Ir -Git) 









-Ir _ ^ L 

2 IP 2r„ 





TV 
2 








■7p 



^Total ) 



(48) 



2tJ 



with the initial condition 



Z^Total(i = 0) 



/ 1 

10 

10 

V 1 / 



(49) 
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The solution of the above differential equation is 



^Total(i) — 



f G{0,t) -fr Jo G~\0,t')Git',t)dt' 
e-'''-*G-\0,t) 

6"^*+'''' + ' 

\ 








-(*+7p+27r-)* 



(50) 



where 

G(ii,t2) =exp|^'G(t)dt|. (51) 

Comparing this matrix with the one in Eq. ()44|) . one can 
see that the full dynamics between times and t is equiv- 
alent to a measurement that is followed by relaxation and 
dephasing operations. With some straightforward alge- 
bra, one can extract the measurement fidelity: 

p ^ !o2Git')e-'^''''G'HO,t')dt' 
2 - /q 2G(i')e-^"*'G-2(0, t')dt' 





G-2(0,t') 




2 + /f""'e-'^-*'d 


G-2(0,tO 



(52) 



Although the above expression looks rather complicated, 
one important result can be seen relatively straightfor- 
wardly. The only difference between this expression and 
the one that would be obtained in the absence of de- 
coherence is the factor e^''"'^* . This factor represents the 
rather intuitive fact that the parts in the QPC signal that 
are recorded at later times carry less value for purposes 
of determining the initial qubit state than those recorded 
at earlier times, with an exponentially decaying function. 
What is less intuitively obvious is how the exponential- 
decay function enters into the expression for the fidelity, 
as can be seen from Eq. ([5^ . 

It is worth making a comment here regarding the 
'value' of the output signal in the absence of decoherence. 
One might be tempted to say that since the fidelity in- 
crease slows down as time goes by (see Fig. 3), the 'value' 
of the QPC's output signal decreases as time goes by in 
that case as well. If that were the case, this decrease 
in the 'value' of the measurement signal would not be 
related to decoherence. However, one can see that this 
is not the case by considering the simple case /3 = in 
the absence of decoherence. In that case the important 
quantity is the time-averaged QPC current throughout 
the measurement duration, without any weight function 
that decreases in time. All parts of the signal are there- 
fore given equal importance when extracting the mea- 
surement result. 

In order to calculate the average measurement fidelity 
in the presence of decoherence, one must average the 
measurement fidelity over all different input states and 
measurement outcomes. Since any calculation involving 
the above expressions seems to be a daunting task, we use 
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Td/ r,n 

FIG. 8: The maximum achievable fidelity as a function of the 
ratio between the decoherence and measurement times. The 
red dots are obtained using numerical simulations, and the 
solid line is given by Eq. (f55)l . 

the following hand-waving calculation. In the absence of 
decoherence, the measurement fidelity after duration t 
can be expressed as 



F{t) w 1-exp 




= / — e-*/""di. (53) 

Jo "^m 

Decoherence reduces the 'value' of the later parts of 
the measurement record with an exponentially decaying 
function (with characteristic time Td). An estimate for 
the average fidelity in the presence of decoherence can 
therefore be given in the form 

ft 1 

F{t) « / — e"*/^'" X e-*/^<^di. (54) 

JO "^m 

Even if the measurement time is taken to infinity, where 
we have _F 1 in the absence of decoherence, we find 
that the fidelity only reaches the value 



Note, in particular, the simple limits that -Fmax when 
Td ^ Tm, and i^max ^ 1 whcu r<j > Tm. We have per- 
formed numerical simulations of the measurement pro- 
cess in the presence of relaxation, and the result are 
shown in Fig. 8. The results agree reasonably well with 
those of the hand-waving argument given above. 
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VI. CONCLUSION 

In conclusion, we have analyzed the question of what 
information can be extracted from the output signal of a 
detector that weakly probes the state of a qubit when the 
qubit Hamiltonian induces oscillations between the two 
eigenstates of the probed operator. We have shown that 
the measurement basis is determined stochastically every 
time the measurement is repeated. In the case of weak 
qubit-detector coupling, the possible measurement bases 
cover all the possible directions. The measurement basis 
and the measurement result can both be extracted from 
the detector's output signal. We have also shown that 
the information acquisition rate is almost independent of 
the angle (3 between the direction defining the probed 
basis and that defining the qubit Hamiltonian. In other 
words, given enough time, the detector will produce a 
high-fidelity measurement result, regardless of the value 
of /3. These results show that, under suitable conditions 
and by proper analysis, the detector's ability to obtain 
high-fidelity information about the state of the qubit is 
not affected by the apparent conflict between the mea- 
surement and coherent-precession dynamics. 

In the course of our analysis, we have introduced an 
equation that supplements the stochastic master equa- 
tion for purposes of extracting measurement information 
from the observed signal. We have also introduced a pro- 
cedure for performing quantum state tomography that is 
well suited for the case where the different measurement 
bases do not have a simple relation to one another. 

We have outlined the generalization of our analy- 
sis to the case where the qubit experiences additional, 
non-measurement-related decoherence. In this case, the 
'value' of the detector's output signal decreases with the 
time at which it is recorded. We have demonstrated this 
fact by analyzing the case where the measurement is per- 
formed in the presence of relaxation. 
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Appendix A: The QPC-current probability 
distribution for a short time interval of size 5t 

In this Appendix, we derive the expression for the 
probability distribution P[/, 51, 5t\ that governs the 
stochastically determined values of the QPC current /. 
We use the discretization parameters 51 and 5t: SI defines 
the size of a finite interval of values for / that we identify 



with a single value of /, and 5t is the size of the time 
interval under consideration. In the derivation below, we 
assume that the qubit is in one of its charge states (for 
definiteness \L)) and that there is no mixing between the 
charge states. We denote the average value of the cur- 
rent by II for the state \L) and Ir for the state \R). 
We are assuming that the probability distribution can 
be approximated by a Gaussian distribution function: 



PL[I,5I,5t]ocexp{-a{l - Il^Y 



(56) 



In order to determine how 5t enters in the expression 
for [I, 51, 5t] we consider the effect of performing N 
measurements in a row. For sufficiently small dt, the 
results of any calculation should not depend on whether 
we treat these measurements as N independent measure- 
ments or as a single measurement of length Ndt. Since 
the values of / obtained in the different measurements 
are independent, their average will have a standard devi- 
ation a/^/N, where a is the standard deviation of / for 
an interval of length 5t. In other words, the standard de- 
viation in / must be proportional to 1/ \/5i. Since the co- 
efficient a in Eq. ([56|) is related to the standard deviation 
cr by (T = l/\/2a, we find that a must be proportional to 
5t: a = a5t. 

We now introduce the measurement time Tm by re- 
quiring that if 5t were to reach r^, the standard de- 
viations (or in other words, the widths) of the proba- 
bility distributions PL[I,5I,5t\ and PR[I,5I,5t\ will be 
equal to half the distance between their center points: 
1 / \/2aTm = \Ir — Il\I'2. Note that this definition does 
not necessarily agree with other definitions of the mea- 
surement time that appear in the literature. Assuming 
for definiteness that Ir> II and defining A/ = Ir — Il, 
we find that 5 = 2/ [t„i [A/j^^ Using the normaliza- 
/ Pl[I , 51 , 5t]dl / (51) = 1, we find 



tion J^i PL[I,5I,5t] 
that 



PL[I,5I,5t] 



25t {5iy 
\\ 7rr„ (A/)' 



■ exp ■ 



25t {I - 



A/ 



(57) 

This expression is used for deriving the matrices 
Um[I , 51, 5t\ in the main text of the paper. 

Appendix B: The relation between a given 
output signal and its 'opposite' 

In this appendix we would like to establish the rela- 
tion between the measurement results of a certain signal 
and its 'opposite'. If a signal is given by Ii{t), with the 
time t running from to i/, the opposite signal is given 
by l2{t) = I — I\{t). In other words, the opposite sig- 
nal is obtained by taking the mirror image of the signal 
about the central current value /. We start by stating 
the relation: if 
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?7Totai[/i(i -.O^tf), 51, 6t] - C/Rot[/i(t : ^ tf),SI, St] x C/Meas[/i(i : ^ t/), 61, 5t], (58) 

with 

C>Moas[/i(i :0^tf),SI,St] = iVi) (^il + /f^|V'2> (^2| , (59) 

then 

UTot^iMt : ^ t/), SI, St] = ;7Rot[/2(t : i/), (ff, x UucMt ■ -> t/), SI, St], (60) 

with 

?/Moas[/i(i :0^tf),SI,St] = IV2) (^2| + (^il , (61) 



and URotMt : ^ i/),a^,(5i] = i7Rot[/i(i : ^ 
tf),SI,St]. 

For = 0, both C/Totai[/i(i : ~* tf),SI,St] and 
C^Totai[^2(i : — > t/), (ff, are given by the unit matrix, 
which obeys the relation given above. We now assume 
that the relation holds for a given value of , and we try 
to establish that it will continue to hold at an infinites- 
imally later time tf + St. The two changes that can oc- 
cur during this infinitesimal interval are a measurement- 
induced weak projection or a Hamiltonian-induccd uni- 



tary transformation. Since Unotihit : tf),SI,St] = 
URot[Ii{t ■ — > tf),SI,St], apphcation of exp{—iHSt} 
to [/Totai[/i(i : ^ tf),SI,St] and UTot^hit : ^ 
tf),SI, St] clearly cannot affect the above relation between 
the total evolution matrices. Demonstrating that appli- 
cation of the weak projection also does not affect the rela- 
tion is less straightforward. Using Eq. (fTU]) and introduc- 
ing the (infinitesimal) quantity q — 5t[Ii[tf )—I]/[Tm^I), 
we find that 



C/Totai[/i(i ■(i^tf+ St), SI, St] cx (1 + qa^) X UKot[h{t : ^ tf),SI, St] x 1^2) (^2! + V^2 IV'i) (^i| j 

= UKot[h{t ■■ ^ tf),SI,St] X (l-fgafi) x (v/i^|^2) m + V^I^Ai) (^i|) ,(62) 



where 



•^ft = ^Rot[^l(^ ^ ^/)' '2"' X (1 + 9'^") X UKot[h{t ■■ ^ t/), <ff, St]. 



(63) 



We now introduce the Pauli matrices = \ — IV'2) (V'2| and a^, for which (/(Jn = ^j^CTj; + (/zCTx. Defining 

ap = (\/Fr+ V7¥)/2 and a„ = {^/Pl ~ ^fK)/2, we find that 



C>Totai[/i(i:0^t/ + ^t),a",,5t] cx UKot[h[t 

= ?/Rot[/i(i 

= UKot[h{t 



^ t/), (5/, St] X (H- -f g^CT^) X (op -f a^o-^) 

i/), (5/, (5i] X (flp + apQ^a-x + a^qzcr^ + a^cTz ~ ia^Qx^y + arnQz) 

t/), (5/, St] X exp{-mmq:rO-y/ap} 



X [op -I- amq^] -I- 



CLpQx 



CTa; -I- [flm -I- flpg^Jo-^ I + O(g^) 



where 



= C>Rot[/i(t : ^ i/ + St), SI, St] X ;7Mcas[/i(i : t/ + (5i), ff, Ji], 
J7Rot[/i(i ■.Q^tf + St), SI, St] = URot[Ii{t : -> i/), (5/, St] x exp{-m,„(7a;a-j^/ap} 

Op = Op + amQz 



(64) 



CLpQx 



am + apQz. 
I — 



(65) 



We can now carry out the same derivation as above for the 'opposite' signal. The only changes are that we re- 
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place q by —q (therefore replacing q^ by —q^ and q^ 
by —qz) and reverse the positions of Pi and P2 (this 
latter change changes Um into —a™ and leaves Up un- 
changed). These changes leave the infinitesimal rota- 
tion angle amqx/o-p unchanged, and we therefore find 
that the relation Unotihit : ~* tf + St), 51, St] — 
URot[Ii{t : tf + St), SI, St] holds. Similarly a'p is 
unaffected by the above changes. The quantities a'„-^ ^ 
and a^^, on the other hand, change their signs. By 
looking at the expressions for the ap and a,„ in terms 
of the Pi and P2 above, we can see that reversing the 
sign of Qm reverses the positions of the symbols Pi 
and P2 in the matrix t/ivicas, which is exactly the dif- 
ference between UMciis[Ii{t : ^ tf + St), SI, St] and 
C^Mcas[^2(^ ■ ^ ^ tf + St), SI, St] that was stated at the 
beginning of this appendix. We have thus established the 
relation between the matrices UTota\[Ii{t '■ tj), SI, St] 
and UTotMt-.O^tf), SI, St]. 

Appendix C: The measurement fidelity for two 
consecutive weak measurements 

In this Appendix, we show that the measurement fi- 



delity for two misaligned weak measurements can be 
higher than that obtained when the measurement axes of 
the two measurements are aligned with each other plj ). 

With no loss of generality, we take the measurement 
axis of the first measurement to be the z axis (i.e., the 
basis {|0) , We assume that the measurement pro- 

duces one of two possible outcomes. The measurement 
matrices can therefore be expressed as 



^ f ^Ri{l + ei) 



^Riil-ei) 



Ul,2 



v/1 - i?i (1 + ei) 



Vi-^i(i-ei) 



(66) 



where we have neglected any coherent component 
in the measurement matrix, as explained in the 
main text. The above matrices satisfy the ba- 
sic requirement that Ul lO 1,1 + t^i 2^1,2 — 1- 
The second measurement is performed in the basis 
{cos(6'/2) |0) + sin{0/2) |1) , sin(6'/2) |0) - cos(6'/2) |1)}. 
The corresponding measurement matrices are given by 



VR2 (1 + £2) 



U2.1 

U2.2 = v/l-^2(l + e2) 



cos 2 sm 2 cos 2 \ , r^—n ^ 

■ e e ■29]+ \^2 (1 - £2) 
sm I cos I sm | 
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cos 2 sm 2 cos 2 
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sm 



2 e 



sm 2 cos 2 
cos^ I 



(67) 



The combined measurement has four possible measurement matrices, each one corresponding to one of the four 
possible outcomes. For example, for the outcome identified with the index 1 in both measurements, the measurement 
matrix is given by V2.\V\,\- In order to obtain the measurement fidelity for this outcome, we need to consider the 
matrix 



V\iUliU2,iVi.i 











V^i(l-ei) 



i?2 + e2COs6' R^il2^ 

R2G2 sin 6 o 
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62 cos 9 











r 



If we denote the eigenvalues of this matrix by Pi and 
P2 (with Pi > P2), the measurement fidelity for this 
outcome is given by (Pi — P2)/{Pi + -Pj)- For purposes 
of calculating the overall fidelity of the measurement se- 
quence, we take the maximally mixed state. 



-'max. mixed 



(69) 



as the initial state. With this initial state, the probability 
that the above outcome is obtained in a given experimen- 
tal run is given by (Pi -I- P2)/2. Thus, the overall fidelity 
(taking into account all four possible outcomes) will be 
given by the sum of four terms, each of which is given by 
(Pi — P2)/2 for one of the possible outcomes. For Eq. 



straightforward algebra gives the result that 

P1-P2 



= P1P2 V(l + £1^2 cos0)2 - (1 - el){l - el). 

(70) 

Similar expressions can be obtained for the three other 
possible outcomes, resulting in the overall fidelity 



F — RiR2\^ 


/(I 


+ 6162 C0S6')2-(1- 


£?)(!- 


-el) + 


R'iR2\^ 


/(I 


-heie2C0s6')2-(l- 


ef){l 


-eD + 


RiR'2\ 


/(I 


+ eie^ cos6')2 - (1 - 


£?)(!- 


-e'i) + 


R'iR'2\ 


Ai 


+ e'ie'2 cos6')2-(l- 


e?){l 


-£?),(71) 


where P^- = 1 — 


Rj 


and e'j — —ejRj/{l 


~R,). 
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The measurement fidelity can liave its maximum at 
any point between zero and tt (with the proper choice of 
the parameters Rj and ej ; we have verified this statement 
numerically). For the special case Ri = R2 = 1/2 and 
ei = e2 = e > 0, the overall measurement fidelity is given 

by 



e2cos6')2-(l-e2)2 + 



V'(l-e2cos6')2-(l-e2)2j . (72) 

Differentiating this expression twice shows that it has a 
maximum at cos 9 = 0. In fact, in the limit of weak 
measurement, i.e. e <C 1, one can see easily that the 
above expression can be approximated by 



cos - 
2 



+ 



sm - 
2 



(73) 



which has minima at ^ = and 9 = n {F = e) and 

a maximum at 6* = 7r/2 (_F = \/2e)- The relative dif- 
ference between the minimum and maximum fidelities is 
approximately 30%. 

The result that the fidelity of two weak measurement 
can be enhanced when they are not aligned with one 
another offers some explanation for the result of Sec. IV 
that the fidelity increases more rapidly when the qubit 
Hamiltonian and the probed operator do not commute. 
The qubit Hamiltonian causes a rotation in the qubit's 
state while the state is being measured along a fixed axis. 
This situation is, in some sense, equivalent to one where 
a number of different qubit operators are sequentially 
probed. 

A final comment that is worth making here regards the 

possible bases of the combined measurement. As can be 
seen from the analysis of this appendix, when the two 
performed measurements correspond to directions in the 
x-z plane, all the resulting matrices will be real. As a 
result, all the possible effective measurements will be in 
the x-z plane as well (note that a projector along the y 
axis would be complex). In other words, when successive 
measurements are made in one plane, they can only re- 
sult in effective measurements made in the same plane. 
This result explains why in the special case /3 = n/2 all 
the possible measurement bases lay in the x-y plane; suc- 
cessive measurements along the x axis separated by state 
precession about the z axis are equivalent to a sequence 
of measurements that are all performed in the x-y plane. 

Appendix D: Quantum state tomography 

In this Appendix, we derive the function that was used 
to perform quantum state tomography in Sec. IV. 

We start with the objective of performing quantum 
state tomography using the available measurement data, 



which we assume has been obtained using a (generally 
weakly coupled) QPC as explained in the main text. We 
note that since the measurements in the different experi- 
mental runs are generally performed in different bases, we 
look for a procedure that treats all the measurements on 
equal footing and does not have a preferred measurement 
basis, or set of bases. One obvious approach to follow is 
to construct a function that assigns a penalty for any de- 
viation by the ' guess' density matrix Pgucss 

from a given 

experimental result (which takes the form of a measured, 
pure state). The guess density matrix is characterized by 

the spherical coordinates rguess, ^guoss and (/"guess, whereas 
the measured state is characterized by only the angles in 
spherical coordinates 9j and (pj (the index j denotes the 
different experimental runs). The function to be min- 
imized for purposes of state tomography will therefore 
have the form 



T(r, 



guess 5 ^guessi ^guess 



3 



guess 5 ^guess: 0guess; ^ji 4^j) 



(74) 

In order to proceed further, we now consider what is 
done in the standard version of state tomography where 
the average values ((Tj;), (a-y) and (ctx) are measured and 
the density matrix p = (1 + {ax) + {py) Oy + ((Tz) 0^)1'^ 
is inferred from the measurement data. If we take the 
measurement data for a^, we know that we want the 
tomography procedure to produce the result rcos^ = 
(7V+1 - 7V_i)/(7V+i + iV-i), where iV+i and N^i are the 
numbers of times that the measurement results and 
— 1 were obtained, respectively. We can therefore see that 
the function /(rgucss, ^'gucss, '/'guess, ^'j = ^Aj) will only 
be a function of (fJ^, guess) = Tguess cos 6'guess (or, in other 
words, the projection of the vector representing the guess 
density matrix along the measurement axis). Taking this 
measurement data, we find that the function T has the 
form: 

7'((6-2,guess) , A^-M,^-l) = Af-M/((o-2,guess)) 

+iV_i/(- (az,guess)).(75) 

In order to minimize the function T, we take its deriva- 
tive with respect to (az.guess): 



,gucss 

d- ((i^z tguess) 



-^+l/'((o'^,guess)) 



If we now use the Taylor expansion 

00 

/ ((<^2, guess)) ~ ^ ^ guess) , {^^^ 



n=0 



we find that 
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dT((a,,gue.s),iV+i,iV-i) 



d iaz. 



guess/ 



= X! - ^^-l)Cn (CTz,guess)" + + Af-l)c„+i ((Tz,guess)"^"^ j 

n=0,2,4,... 

= (A^+l+A^-l) ^ {c„ ((T^^corrcct) + C„+i (o-^^guoss)} (^z.g 



:, guess/ J 



n=0,2,4,. 



(78) 



r 



where (o-^.corrcct) = - Ar-i)/(iV+i +iV-i). In or- 

der to ensure that the function T has a minimum at 
(<5-z,gucss) = {cTz, correct), we must choose c„ = -c„+i for 
every even number n. The simplest choice is cq = — ci = 
— 1/2 and C2 = C3 = ... = 0. This choice gives the func- 
tion 



TO 



gUCSSi ^gUCSSi fgucss. 



2:, guess/ 



(79) 



which, when generahzed to the case where the measure- 
ment axis is not necessarily the z axis, gives Eq. ([57|) . 



The detailed analysis of the different possible choices 
of c„ in the function T, and the stability and reliability of 
these different choices, is beyond the scope of this paper. 
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